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Abstract 

In this paper, we calculate the decay rate and CP asymmetry of the B s — > ir + ir~ decay in 
perturbative QCD approach with Sudakov resummation. Since none of the quarks in final states 
is the same as those of the initial B s meson, this decay can occur only via annihilation diagrams 
in the standard model. Besides the current-current operators, the contributions from the QCD 
and electroweak penguin operators are also taken into account. We find that (a) the branching 
ratio is about 4 x 10 -7 ; (b) the penguin diagrams dominate the total contribution; and (c) the 
direct CP asymmetry is small in size: no more than 3%; but the mixing-induced CP asymmetry 
can be as large as ten percent testable in the near future LHC-b experiments. 

PACS numbers: 13.25.Hw,12.38.Bx 
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I. INTRODUCTION 



In recent years, rare B decays are attracting more and more attentions, since they 
provide a good opportunity for testing the Standard Model(SM), probing CP violation 
and searching for possible new physics beyond the SM. Since 1999, the data sample of 
the pair production and decays of B mesons collected by BaBar and Belle Collaborations 
is increased rapidly In the future LHC-b experiments, B s and Be mesons can also be 
produced, and the rare B decays with a branching ratio around 10 -7 can be observed. 
The rapid progress in current B factory experiments and the bright expectation in LHC-b 
experiments induced a great interest in the studies of rare decays of B meson. 

The rare decay B s — > 7r + 7r~ can occur only via annihilation diagrams in SM because 
none of quarks in final states is the same as those of the initial B s meson. The usual 

n 

method to treat non-leptonic decays of B meson is Factorization ApproachfFA) which 
has achieved great success in explaining many decay branching ratios 0, 0, Q| . However, 
this method failed in describing B s — > tt + tt~ decay, because we need the 7r — > 7r form 
factor at very large momentum transfer O(Mb)- So far, little is known about the form 
factor at such a large momentum transfer in FA. In the QCD factorization approach j^, 
one cannot perform a real calculation of the annihilation diagrams, but estimating the 
annihilation amplitude by introducing a phenomenological parameter. In this paper, we 
calculate the branching ratio and CP asymmetries of B s — > 7t + tt~ decay by employing 
the perturbative QCD approachfPQCD) p. This method has been developed for the 
studies of the B meson decays |7[ and successfully applied to calculate the annihilation 
diagrams 0, |^. When the final states are light mesons such as pions, the perturbative 
QCD approach(PQCD) can be safely used because of asymptotic freedom of QCD jlfl| |. 

In the next section, we give our theoretical formulas for the decay B s — > 7r + 7r~ in 
PQCD framework. In section 3, we give the numerical results of the branching ratio of 
B s — > 7r + 7r~ and discuss CP asymmetry of the decay. 

II. PERTURBATIVE CALCULATIONS 



The related effective Hamiltonian for the process B s — > 7r + 7r is given by |9j, 

h cS = -4{v: h v us [Ci^Ox^) + cmom] - v; b v ts £ c^o^)}, (i; 
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where Cj(//)(i = 1, • • • ,10) are Wilson coefficients at the renormalization scale /x and the 
operators Oi(i — 1, • • • , 10) are 

Oi = (SiUj)v-A(Ujbi) V -A, 

3 = {sibi) V - A ^2{qjqj)v-A, 

q 

4 = (sibj)y-A ^(qjq^v-A, 

i 

5 = (sibi) v ^A^2(qjqj)v+A, 

q 

6 = (Sjbj) V - A y^(gjgj)y+A, 

9 

3 

07 = ^(sibj)v-A ^ e g (g j g j ) y+A , 

3 

°s = -(s i &j)y-A^e (? (g j gj)y + A, 
3 

9 

3 

O10 = ^ibjW-A^eqfaqdv-A- (2) 
9 

Here i and j are S'C/ (3) color indices, the sum over q runs over the quark field that are active 
at the scale \i = 0(m&), i.e., q G {u, d, s,c, b}. Operators 0\,Oi come from tree level, 
3 , 4 , 5 , 6 are QCD-Penguins operators and 7 , 8 , O9, Oi come from electroweak- 
Penguins. 

In the PQCD approach, the decay amplitude is separated into soft (<£>), hard(if), and 
harder(C) dynamics characterized by different scales. It is conceptually written as the 
following, 

Amplitude- Jd%d 4 k 2 d% ^[Cit^Bih^ih^AhWkuh,^)], (3) 

where fcj's are momenta of light quarks included in each mesons, and Tr denotes the trace 
over Dirac and color indices. C(t) is Wilson coefficient which results from the radiative 
corrections at short distance. $m is wave function which describes the hadronization of 
mesons. The wave functions should be universal and channel independent, we can use 
$m which is determined by other ways. The hard part H is rather process-dependent. In 
the following, we start to compute the decay amplitude of B s — > it + it~ . 
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Since we set B s at rest, in the light-cone coordinates, the momentum of the B s , ir 
and 7r + are written as : 

Mr Mr Mr 

Pi = ^§(1, 1, Or), P 2 = 0, Or), P 3 = ^§(0, 1, Or). (4) 

Denoting the light (anti-)quark momenta in B, tt~ and tt + as kx, k 2 and k 3 , respectively, 
we can choose: 

h = (xipl,0,k 1T ), k 2 = (x 2 p 2 ,0,k 2T ), h = (0,x 3 p^ ,k 3T ). (5) 

According to effective Hamiltonian(|T|) . we draw the lowest order diagrams of B s — ► 
7r + 7r _ in FigO For the factorizable diagrams (a) and (b), we find their contributions 
cancel each other, which is a result of exact isospin symmetry. For the non-factorizable 
diagrams (c) and (d), the contribution comes from tree operator is 

j rl />oo 

Mj = QAtcCfMb / dx\dx 2 dx 3 / b\db\ b 2 db 2 4>B(xi,b\) 

C 2 {ti)a s {t 1 ){-x 3 ^{x 2 )^{x 3 ) - r\x 2 + x 3 )^{x 2 )<f)^{x 3 ) 

+ r 2 (x 2 - x 3 )(j)P{x 2 )(j)l{x 3 ) + r 2 (x 2 - x 3 )(p^(x 3 )(j)l(x 2 ) 
- r 2 (x 2 + x 3 )(pl{x 2 )(t)l{x 3 )}h^\x l ,x 2) x 3l b l ,b 2 )E m (ti) 
+ C 2 (t 2 )a s (t 2 ){x 2 (j)^(x 2 )(f)^(x 3 ) + r 2 (2 + x 2 + x 3 )^(x 2 )^(x 3 ) 
+ r 2 (x 2 - x 3 )4% (x 2 )(frl(x 3 ) + r 2 (x 2 - x 3 )(j)l(x 2 )(j)^(x 3 ) 

+ r 2 (-2 + x 2 + x 3 )4>l{x 2 )4>l{x 3 )}h ( i\x 1 , x 2 , x 3 , b u b 2 )E m {t 2 )] , (6) 

where r = m 0n /mB = m^/[mB(Tn u + ma)]- Cf = 4/3 is the group factor of the 577(3) c 
gauge group. The expressions of the meson distribution amplitudes 0m, the Sudakov 
factor E m (ti) and the functions ha ' (xi,x 2 ,x 3 ,bi,b 2 ) are given in the appendix. The 
contribution of penguin- diagrams M^ can be obtained by replacing C 2 (ti)(i = 1, 2) with 

a p {U) = 2C 4 (*i) + 2C 6 (U) + ~C 8 {U) + ~C 10 {ti), (7) 

in Eq.(jBJ). The explicit expressions of QCD corrected Wilson coefficients C 2 , C4, C$, C$ 
and C10 as a function of scale t can be found in the Appendix of Ref . . 
Now, the total decay amplitude for B s — > 7r + 7r~ is given by 

A = V: b V us Ml - V t lV ts M? = V: h V us M T a [l + ze^\ (8) 
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FIG. 1: The lowest order diagrams for B® — ► -k + tt decay. 



where z = \ V t * b Vt s /V* b V us \ \M£ /Mj\, V ub ~ \V U b\e %1 and 5 is the relative strong phase 
between tree diagrams(Mj) and penguin diagrams(M ( f ). z and 5 can be calculated from 
PQCD. 

The decay width is expressed as 



r(3 



7T 7T 



G%M B lA]2 = ^k KbVuaM T ?[1 + Z 2 + 2zcos{6 _ 7)] _ (9) 



1287T 1 ' 1287T 

Similarly, we can get the decay width for B° s — > tt + it 



where 



A = V ub V: s Mj - V tb V t *M = V ub V* s M'+'{l + ze 



IT IT 



128vr 1 1 ' 



(10) 



11) 



III. NUMERICAL EVALUATION AND SUMMARY 



The following parameters have been used in our numerical calculation: 

M Bs = 5.37 GeV, m 07r = 1.4 GeV, A£= 4 D = 0.25 GeV, f Bs = 236 MeV, 



U = 130 MeV, t b o = 1.46 x 10~ 12 s, \V t * b V ts \ = 0.0395, \V* b V, 



0.0008. 



(12) 



5 



We leave the CKM phase angle 7 as a free parameter to explore the branching ratio and 
CP asymmetry parameter dependence on it. In SM, the CKM phase angle is the origin 
of CP violation. 

From Eq.(jnj) and we get the averaged decay width for B® (B®) — ► 7r + 7r~ 

r 2 A/f 3 

r(B° s m - tt+O = %^(|A| 2 /2 + |A| 2 /2) 

12o7T 

= -^\V: b V us Mj\ 2 [l + 2zcos 1C os6 + z\ (13) 

Using the above parameters, we get z = 13.4 and 5 = 168° in PQCD. Eq. fT3|) is a 
function of CKM angle 7. In Fig. we plot the averaged branching ratio of the decay 
B®(B®) — > 7r + 7r~ with respect to the parameter 7. From Fig. EJ we can see that [lsj ]: 

Br(5°(£°) -> 7T+7T-) = (4.2 ± 0.6) x 1(T 7 , (14) 

for < 7 < 77. The number z = \V*V ts /V* b V us \ |Mf/Mj| = 13.4 means the amplitude of 
penguin diagrams is about 13.4 times more than that of tree diagrams. Therefore almost 
all the contribution comes from penguin diagrams in this decay and the branching ratio 
is not sensitive to 7. 

In Ref . . Beneke et al have estimated the branching ratio for B s — > 7r + 7r~ in QCD 
Factorization approach. In order to avoid the endpoint singularities, they introduced 
parameters to replace the divergent integral. In this approach, they estimated that the 
branching ratio of this decay is (0.24— 1.55) x 10~ 7 with those phenomenologicalparam- 
eters. In our work, the calculation has no endpoint singularity because of fcrj6j. Our 
predicted result is larger than their simple estimation, which can be tested by the exper- 
iments. 

For the experimental side, we notice that there is only upper limit of the decay B® — > 
7t + tt~ given at 90% confidence level [3] 

Br(B° s -> tt + tt-) < 1.4 x 10~ 4 . (15) 

Obviously, our predicted result is still far from this upper limit. 

In SM, CP violation comes from interference between amplitudes with different CP 
eigenvalues. The strong interaction eigenstates B® and B® can mix through weak inter- 
action, i.e. B®-B® oscillation. By experimental observation we can know whether CP is 

nn 

conserved. For the B®-B® system, the CP asymmetry is time dependent |3l. Il3| : 

A CP {t) ~ Ac r P cos(Amt) + a e+e > sin(Amt), (16) 
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FIG. 2: The CP-averaged branching ratio of B®(B S ) — > -k + it decay as a function of CKM angle 
7- 

where Am is the mass difference of the two mass eigenstates of mesons. Af'p is the 
direct CP violation parameter while a e+e > is the mixing-related CP violation parameter. 
The direct CP violation parameter is defined as 

Adir = T(B° S -» 7T+7T-) - r(g -> 7T+7T-) = 2z sin 7 sin j 

CP r(Sj -»• 7T+7T-) + T{B° -»• 7T+7T-) 1 + 2z COS 7 COS 5 + Z 2 ' 1 ' 

Using Eq.@ and (jlOj). we can compute the parameter The direct CP asymmetry 

Aq P has a strong dependence on the CKM angle, as can be seen easily from Eq.(|17|) and 
Fig. El From this figure one can see that when the CKM angle 7 is around 7r/2 the direct 
CP asymmetry reaches its peak, which is about 3%. The small direct CP asymmetry is 
also a result of small tree level contribution. 

The mixing-related CP violation parameter in Eq. f)lf)j1 is defined as [^] 



— 2Im(Xcp 

a e+e' = 

where 



<w - TT^p US) 



^-v^mtQWY (19) 

In Fig. EJ we study the mixing CP violation parameter a 6+6 > of the decay B s — > 7t + 7t~ as a 
function of CKM angle 7, just like the case of direct CP violation, it is almost symmetric 
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FIG. 3: Direct CP violation parameters of B s — > ir + ir decay as a function of CKM angle 7. 

and the symmetry axis is near 7 = tc/2. its peak is close to —14.5%. The possible large 
CP asymmetry might be observed at LHCb experiment in the future, this would help us 
to determine the value of CKM angle 7. 

In conclusion, we study the branching ratio and CP asymmetry of the decay B° s (B° s ) — > 
7r + /T~ in PQCD, we find that the branching ratio is at the order of 10~ 7 and there are large 
CP asymmetries in the process, which may be measured in the future LHC-b experiments 
and BTeV experiment at Fermilab. This small branching ratio, predicted in the SM, make 
it sensitive to possible new physics contribution. 
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FIG. 4: Mixing-related CP violation parameter a e+e ' of B®(B S ) — ► 7r + 7r decay as a function of 
CKM angle 7. 

APPENDIX A: SOME FORMULAS USED IN THE TEXT 

For 5, meson wave function, we use the same wave function as B meson , despite 

the possible SU(3) breaking effect 



The parameter Uf, = 0.4 GeV, is constrained by other charmless B decays |9j, |14|. And 
Ns a = 114.0GeV is the normalization constant using fs a = 236MeV. 

The 7T meson's distribution amplitudes are given by light cone QCD sum rules 1 : 



(fi Bs (x,b) = N Bs x 2 (l -x) 2 exp 



" M 2 B x 2 



(Al) 




(A2) 



where t — 1 — 2x. The Gegenbauer polynomials are defined by: 



Cl'\t) = i(3t 2 - 1), C\' 2 {t) = i(35t 4 - 30t 2 + 3), 
Cl'\t) = 3 -(5t 2 - 1), Cl'\t) = ^(21t 4 - Ut 2 + 1). 



(A3) 



Since the hard part is calculated only to leading order of a s , we use the one loop 
expression for the strong running coupling constant in our numerical analysis, 

47T 



(A4) 



where (3q = (33 — 2nf)/3 and nf is the number of active quark flavor at the appropriate 
scale /i. A is the QCD scale, we set A = 250MeV at rif = 4. 
The function E m (ti) in Eq.© are defined by 



Efn, (ti 



-S B (ti)-S w+ {U)-S^{U) 



(A5) 



where S B , S w +, S^- result from summing both double logarithms due to infrared gluon 
corrections and single ones caused by the renormalization of ultra-violet divergence. They 
are defined as: 

" dft 



S B (t)=s(x 1 P+,b 1 ) + 2 



-7(a s (/x)), 



S n -{t) = s(x 2 P 2 + ,b 2 ) + s((l-x 2 )P 2 + ,b 2 )+2 f ^(a s (fl)) 

S w+ (t) = s(x 3 P 3 -,b 3 ) + s((l-x 3 )P 3 ,b 3 ) + 2 [ %(a s (jl)), 

where s(Q,b) called Sudakov factor is given as j3] 

Q 



s(Q,b) 



t/t 



In 



with 



L>F — 




7T 




2a s 


In 


3 TT 





67 7T 
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A(a(/j)) + B(a s (fi)) 



9 

D 2 7 E-1 



fit H — /3n In 
3 27 f 3 f 



-DE 



(A6) 
(A7) 
(A8) 

(A9) 

(AlO) 
(All) 



here 7^ is the Euler constant, nf is the active flavor number. 

The functions h^a\ and in Eq.© come from the Fourier transformation of propa- 
gators of virtual quark and gluon. They are defined by 

h ( a J \xi,x 2 ,x 3 ,b 1 ,b 2 ) = 



^R i o\M By ^x^b 1 )3 (M B ^x^ 3 'b 2 )e(b 1 - b 2 ) 



+ (bi <-> b 2 



KoiMsF^h), for Fg> > 
m^\M Bx /\F^\b 1 ), for F 2 <0 



(A12) 
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where H. 



(z) = Jq(z) + iYo(z), and F^s are defined by 



Ffa = xix 3 - x 2 x 3 , F ( 2 2) = xi + x 2 + x 3 - xix 3 - x 2 x 3 . 



(A13) 



The hard scale t{(i = 1,2) in Eq.© are taken as the largest energy scale in the H to 
kill the large logarithmic radiative corrections: 
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